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Abstract Generally, the first step in modeling molecular magnets involves obtaining the low- 
lying eigenstates of a Heisenberg exchange Hamiltonian which conserves total spin and belongs 
usually to a non-Abelian point group. In quantum chemistry, it has been a long standing problem 
to target a state which has definite total spin and also belongs to a definite irreducible represen- 
tation of the point group. Many attempts have been made over years, but unfortunately these 



O , have not resulted in methods that are easy to implement, or even applicable to all point groups. 

^ Here we present a general technique which is a hybrid method based on Valence Bond basis and 

the basis of z-component of the total spin, which is applicable to all types of point groups and 
is easy to implement on computer. We illustrate the power of the method by applying it to the 



molecular magnetic system, Cu^Fe^, with cubic symmetry. We emphasize that our method is 
applicable to spin clusters with arbitrary site spins and is easily extended to fermionic systems. 
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2 1 Introduction 
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The field of molecular magnetism has witnessed an explosion in the number of systems that ex- 
hibit molecular magnetic phenomena such as quantum resonant tunneling and photomagnetism 
(see reviews HI |2[). This explosive growth has also presented challenges to theorists model- 
ing these systems. The problems encountered by theorists begin with modeling the nature of 
exchange interactions between pairs of magnetic ions. While an electronic many -body Hamil- 
tonian has to be solved for determining the nature of exchange, this is often circumvented by 
guessing the nature of exchange based on the knowledge of the ligands, electron configuration 
of the transition metal ion and the geometry of the complex. The second problem concerns with 
obtaining the eigenstates of the exchange Hamiltonian, 
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describing the coupling between pairs of magnetic ions (ij) with exchange constant //,-. Of- 
ten, the Fock space of the Hamiltonian of the magnetic system could be very large (in case of 
MnyiAc, it is as large as a hundred million 01 and in Fe\2 ferric wheel, it is more than two bil- 
lion [3J) and obtaining even a few low-lying states of the Hamiltonian could pose a challenge. 
Since the exchange Hamiltonian conserves both total spin and z-component of total spin (Ms), 
the problem can be simplified by specializing the basis, in which the matrix representation of 
the Hamiltonian is sought, to the case of fixed z-component of the total spin. Further simplifi- 
cation could come from exploiting spatial symmetries of the model. An ideal situation would 
correspond to one in which all the spin and spatial symmetries are utilized to construct a fully 
symmetrized basis to minimize the size of the Hamiltonian matrix that needs to be diagonalized. 

The conservation of the Sf ot , the total z-component of spin is trivially achieved by choos- 
ing from the Fock space, states whose total Ms value corresponds to the desired value. This is 
possible since individual S^ operators commute with the S z tot operator. It is also quite straight- 
forward to set up the Hamiltonian matrix in this basis and solve for a few low-lying states in 
cases where the Hilbert space is spanned by a few hundred million states 11301 . However, the 
exchange Hamiltonian of molecular magnets often consists of interactions that are geometri- 
cally frustrated. In such a system, the ground state spin is often not predictable and one needs 
to obtain the lowest state in each total spin subspace to fix the spin of the ground state. Besides, 
low energy states with different total spins lie close in energy and it is numerically difficult to 
achieve convergence to nearly degenerate eigenstates unless they can be dispersed into orthogo- 
nal Hilbert spaces. We can partially alleviate this problem by employing the parity symmetry of 
the exchange Hamiltonian. This symmetry corresponds to rotation of all the spins in the system 
around the x or y-axis by an angle n which leaves the Hamiltonian invariant, in the Ms = 
sector. 



[H,R y (n)\ = 0; where, R y (B) = e~ iQS y/ h (2) 

The action of parity operator (P) on a basis state with site m s values m\,m,2,m?, .... m n , is to flip 
all the spins in the system, i.e., 



P\m\m2. . .m n ) = (— l)^\ — m{ — ni2 . . . — m n ); where, r\ = ^si (3) 

i 

Thus, the parity operator which conserves the total Ms value, only when Ms = 0, is a sym- 
metry element of the Hamiltonian matrix in the Ms = sector. In the general case where the 
individual objects have spin su if Y<i s i ls even men symmetric (antisymmetric) combination of 
the basis states, under parity, will span a space of even (odd) total spin states. The method 
can be extended easily to systems which allow only half odd-integer total spin. Since in most 
cases, the lowest excited state usually has a spin which is one different from that of the ground 
state, this symmetry renders it easy to obtain the spin gap accurately. However, the size of 
the Hilbert space is only reduced by approximately half of the size of the full Ms = space, 



by using this symmetry. Use of parity is still advantageous as exploiting spatial symmetries is 
straightforward, even when the point group is not Abelian. 

Construction of spin adapted configuration state functions (CSF) has been a problem of long 



standing interest in quantum chemistry. The CSFs are simultaneous eigenstates of total Sf ot and 
Sf ot and setting up the Hamiltonian matrix in this basis leads to matrices of smaller size besides 
allowing automatic labeling of the states by the total spin. Besides, the eigenvalue spectrum 
is enriched, since we can obtain several low-lying states in each total spin sector. This is in 
contrast to obtaining several low-lying states in a given total M5 sector as the latter would have 
states with total spin S to t >Ms. There are many ways of achieving this task [5]. The simplest 
method involves setting up the matrix of the total spin operator, Sf ot , in the CSF basis of fixed 
Ms and obtain the eigenstates corresponding to given total spin value; these eigenstates which 
are linear combinations of the constant Ms CSFs are then the spin adapted CSFs. For large 
systems this method is not practical. Another method which is some times used is the Lowdin 
projection method BBEl in which a projection operator P$, given by, 



P s = n s ^s(S 2 -S / (S / +l)) (4) 

is used to project out all undesired states of spins S' except the spin S of interest, from a given 
CSF. The methods that have been extensively in vogue for construction of the spin adapted CSFs 
are the Graphical Unitary Group Approach (GUGA) (HllKTOl, Symmetry Group Graphical Ap- 
proach (SGGA) [11] and the Valence Bond (VB) approach EUnHHEHn). In tne GUGA 
method, the total spin adapted CSFs are represented as Shavitt graphs or Paldus arrays and the 
matrix element of a term (which corresponds to a generator of the unitary group) in the Hamil- 
tonian between any two CSFs is obtained by comparing the two arrays or graphs corresponding 
to the CSFs. Similarly, graphs are used to represent spin adapted CSFs in the SGGA method 
and rules for computing matrix elements between two CSFs of a term in the Hamiltonian have 
been derived (see for example reviews [[T8l[T9ll20l ). Using these methods it is possible to carry 
out large scale configuration interaction (CI) calculations. While in all these methods the total 
spin symmetry is fully exploited, spatial symmetry adaptation is not an easy task [fT8l . The 
CSFs are each built up of several orbitals with each orbital in general transforming according to 
some specific irreducible representation of the point group of the system. The direct product of 
the irreducible representations of a general symmetry group is not a single irreducible represen- 
tation of the same group. Thus it is not possible to associate an irreducible representation with 
a given CSF unless the point group to which the system belongs is an Abelian group ll2Tll22ll . 
Otherwise, symmetry operation on a CSF leads to a linear combination of many CSFs which 
is in general difficult to construct. For small dimensionalities of the Hilbert spaces, matrix rep- 
resentation of the symmetry operators can be obtained in the space of CSFs. The projection 
matrix for a given irreducible representation can be constructed from these matrices and from 
these, the symmetry adapted CSFs [|23l . However, this approach is of limited value in real large 
scale problems lfT8l . In quantum chemical literature this difficulty is bypassed by dealing with 
Abelian subgroups of the systems point group [24] [251 . But, this can lead to ambiguities in 
assigning the irreducible representation of a state [|26l . Among the methods for constructing 



spin adapted CSFs, the VB method is the simplest and will be considered in the next section. 

The ultimate goal of symmetry adaptation is to be able to exploit the full spatial and spin 
symmetries of the system, both for computational efficiency and for labeling of the state by the 
irreducible representation to which it belongs. In the next section, we give a brief introduction 
to the symmetrized valence bond (VB) approach that was developed earlier and highlight the 
difficulties associated with the technique lfT6ll . In the third section, we present a hybrid VB- 
constant Ms method which overcomes these difficulties. In the fourth section, we illustrate an 
application of this method to a magnetic spin cluster. In the final section, we summarize and 
discuss the future prospects for the technique. 



2 The Symmetrized VB Approach 

Exploiting the invariance of both total spin and its z-component is nontrivial, since eigenstates 
of the Sf ot operator expressed as a product of the eigenstates of all the Sj operators are not 
simultaneously eigenstates of the Sf ot operator. The situation is further complicated by the fact 
that in a molecular magnet, often the spins of all the constituent magnetic centers, si are not 
the same. In such a situation, the easiest way of constructing the spin adapted functions is the 
diagrammatic valence bond (VB) method based on modified Rumer-Pauling rules lfT2l[T4ll . In 
this method, a magnetic site with a given spin "s" is replaced by 1st spin-half objects. To obtain 
a state with total spin 5 from N such spin- 1/2 objects from all the magnetic centers, (N — 25) of 
these spin- 1/2 objects are singlet spin paired explicitly, subject to the following restrictions: (1) 
there should be no singlet pairing of any two spin-half objects belonging to the same magnetic 
center (this ensures that the 2si objects are in a totally symmetric combination 11291 ) , (2) a total 
of 25 spin-half objects are left unpaired, (3) when all the spin-half objects are arranged at the 
vertices of a regular polygon with number of vertices equal to number of spin-half objects, N, 
and straight lines are drawn between spin paired vertices, there should be no intersecting lines 
in the resulting diagram and (4) when all the spin-half objects are arranged on a straight line and 
lines are drawn between spin paired objects, these lines should not enclose any unpaired spin- 
1/2 object. These rules follow from the generalization of the Rumer-Pauling rules to objects 
with spin greater than 1/2 and total spin greater than zero. The set of diagrams which obey 
these rules would hence forth be called "legal" VB diagrams. Some legal VB diagrams are 
shown in Fig. 1. A line in the VB diagram between two spin -1/2 objects i and j corresponds to 
the state (oc/(3 ; - — p\a,)/v^2 where we choose a to correspond to | f ) and |3 to | |) orientations of 
the spin. The phase convention assumed is that the ordinal number "z" is less than the ordinal 
number "/'. The 25 spin- 1/2 objects k\ k% % . . . . k2s which are left unpaired can be taken 
to represent the state with Ms = 5 given by a y t 1 OC£ 2 a / t 3 ...CX£ 2s . VB states corresponding to other 
M5 value for this state with spin S, can be obtained by operating the S^ operator on the state 
by the required number of times. Since the exchange Hamiltonian is isotropic, each eigenstate 
in the spin S is (2S+1) fold degenerate and it is sufficient to work in subspace of chosen M5 
value. The VB state corresponding to a given diagram is a product of the states representing 
the constituent parts of the diagram. On a computer, a "legal" VB diagram of any spin can be 
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Figure 1 : Above VB diagram shows spin pairings to yield a total spin S tot =0 state from ten 
spin- 1/2 objects, constituent elementary spins of two spin 1 and two spin 3/2. Its bit represen- 
tation corresponds to unique integer I = 856. The bottom VB diagram shows a S to t=^ state, the 
corresponding unique integer is I = 888. 



uniquely represented by an integer of 2N bits with a one bit representing the beginning of a 
singlet line and a zero bit the ending of singlet line. The unpaired spins are also represented as 
one-bits. Fig. 1 also shows bit representation of typical VB diagrams. 

To spatially symmetrize a VB basis, it is necessary to know the result of a symmetry op- 
erator operating on a legal VB diagram. In general, the resultant of such an operation on a 
"legal" VB diagram is an "illegal" VB diagram. An example of this is shown in Fig. 2 lfT6l . An 
illegal VB diagram can be decomposed into a linear combination of legal VB diagrams, a pro- 
cess that is computationally demanding. In practice, the VB space is broken down into smaller 
invariant subspaces of the symmetry group and within each invariant space, a symmetrized lin- 
ear combination of the VB basis is constructed. However, the structure of the invariant spaces 
is very complex and constructing disjoint invariant spaces is not simple. While constructing 
symmetrized VB basis, a projection operator for a given symmetry representation is employed 
to project the symmetrized linear combinations by acting on each of the VB states in the in- 
variant space. While the number of linearly independent symmetry combinations of a given 
representation is known a priori, the actual linear combinations are obtained by carrying out 
Gram-Schmidt orthonormalization of the projected states. However, since the VB diagrams 
are not orthogonal the orthonormalization process is both computationally involved and time 
consuming. 

Furthermore, in case of molecular magnets containing magnetic ions with spin greater than 
half, the exchange operator between such high-spin centers also generates "illegal" VB dia- 
grams as it involves non-nearest neighbor exchange interactions between constituent elemen- 
tary spins [fT4ll . To illustrate, exchange between a center A with say spin one and a center B 
with spin 3/2, Sa ■ Sb, is expressed as (sa 1 +*a 2 ) ' (*Bi +sb 2 + s b 3 )- These exchange terms op- 
erate on a VB diagram with constituent elementary spins which are non-nearest neighbors and 
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Figure 2: The effect of operation by the C3 symmetry operator about the (1,8) axis. Top left 
shows the initial and final VB diagrams with spin couplings between vertices of the cube shown 
as dark lines. Bottom left shows the same states as spin couplings between vertices of a reg- 
ular octagon. The resultant is an "illegal" diagram. Decomposing the resultant to "legal" VB 
diagrams yields a sum of five VB diagrams shown on the right, with spin couplings between 
vertices on a cube. 



in general generate "illegal" VB diagrams as per the VB rules. Decomposition of the resultant 
"illegal" VB diagrams to "legal" VB diagrams could present a serious bottle-neck in comput- 
ing the Hamiltonian matrix elements. In view of these difficulties, a fully symmetrized VB 
approach to solving Heisenberg exchange Hamiltonian particularly in the context of molecular 
magnets is not feasible. 



3 Hybrid method based on VB Basis and Constant M 5 Basis 

In the constant Ms basis, a basis state of an ensemble of spins s\, S2, • • •, Sff, is represented by 
a direct product of the m s states of each spin such that the total Ms = Y* m i- The basis states 
in the constant Ms representation are orthonormal by construction. Given the definition of a 
line in the VB diagram, every VB diagram can be broken up into a linear combination of the 
constant Ms basis states. A VB diagram with p singlet lines give rise to 2 P basis states in the 
constant Ms basis. To effect the conversion of VB diagrams to constant Ms functions, each 
singlet line gives two states; in one state, the site at which a singlet line begins is replaced by 
an a spin and the one at which it ends by a |3 spin with phase +1 and in the other the spins are 
reversed and the associated phase is -1. Once the VB diagrams are converted to constant Ms 
basis states with constituent spins, it is possible to associate a m ; value with a composite spin, 
given by rm = (n^ — rc/j.)/2, where n^ is the number of up-spin halves and n^ is the number 
of down-spin halves at the site i. However, there is a normalization constant w ; , which follows 
from Clebsch-Gordan coefficients, given by, 



W; 



(2*)! 



1-1/2 



(si + mi)\(si-mi)\ 
6 



(5) 



for a site with spin s, in state m, [|29l . We can assume without loss of generality that the Ms value 
of the VB diagram is also S. Computationally, finding the transformation of a state in the VB 
basis to constant Ms basis is straightforward. We initialize the coefficients in the constant Ms 
basis to zero. We then decompose, sequentially, each VB diagram into constant Ms states and 
update the coefficient of the basis state of corresponding Ms by adding to it the VB coefficient 
times the product of Clebsch-Gordan factors with appropriate phases. The matrix relating the 
VB basis states to constant Ms basis states, C, is a V x M matrix, where V is the dimensionality 
of the VB space and M that of the constant Ms space. 

The construction of the projection operator for projecting all the states of a given symmetry 
representation in a given spin space can now be accomplished. We construct the matrix repre- 
sentation of a symmetry operator, R, of the point group in the chosen spin space by operating 
with R on each state in the constant Ms basis and searching for the resulting state in the list of 
Ms basis states. Each basis state in this representation is carried over to another basis state by a 
symmetry operation of the point group. Thus, the matrix R^ though anMxAf matrix contains 
only one nonzero element in each row; this makes manipulations with this matrix computa- 
tionally fast. The knowledge of the C and the R^ matrices give the effect of operating by the 
symmetry operator R on a VB state as a linear combination of the constant Ms basis states via 
the matrix B#=CRm- The projection operator for projecting out the basis states on to a chosen 
irreducible representation of the point group T is given by, 



Pr = Y,Xr r (R)R (6) 



where, %f r (R) is the character under the symmetry operation R in the character table of the 
point group of the system OTTl . The matrix representation of Pp m the mixed VB and constant 
Ms basis is given by, 



Qr = IXr r (#)B* (7) 

R 

where, Qr is a V x M matrix. However, the rows of the matrix Qr are not linearly indepen- 
dent, since the symmetrized basis spans a much smaller dimensional Hilbert space. The exact 
dimension of the Hilbert space spanned by the system in the irreducible representation T can be 
known a priori. The dimension of the space T, Vr, is given by, 



V r = (d r /h)Y<%(R)%F(R) (8) 



where dr is the dimensionality of the irreducible representation T, h is the number of elements 
in the point group and %(R) is the reducible character for the operation R. The Vr x M projection 
matrix, Pr of rank Vr is obtained by Gram-Schmidt orthonormalization of the rows of the matrix 



Qr until Vr orthonormal rows are obtained. These Vr orthonormal rows represent the projection 
matrix Pp. The M x M Hamiltonian matrix, H^ is constructed in the constant Ms basis which 
is described elsewhere 11301 . Since the basis states in this representation are orthonormal, we 
do not encounter the problem of "illegal" states as with the VB representation. The Vr x Vr 
Hamiltonian matrix in the fully symmetrized basis is given by Pr^M (Pr) an d one could use 
any of the well known full diagonalization routines to obtain the full eigenspectrum or use 
Davidson algorithm to get a few low-lying states of the symmetrized block Hamiltonian in the 
chosen spin and symmetry subspace. 

The above procedure does not lead to the smallest block of the Hamiltonian matrix, when 
the irreducible representation for the block is degenerate, such as the E, T or H representa- 
tions. In such cases, it is advantageous to work with bases that transform according to one of 
the components of the irreducible representation. This can be achieved by choosing an axis of 
quantization and projecting out bases states of the irreducible representation which are diag- 
onal about a rotation about the quantization axes. For example, in the case of the irreducible 
representation that transforms as T, we can choose one of the C3 axes as a quantization axis 
and project the basis states of the irreducible representation using (I+Cg+C^ 1 ) as the projection 
operator. This operator projects states that transform as the Yj component of the three fold de- 
generate irreducible tensor operator. Similarly, for the E representation, we could use any of the 
C2 axis as a quantization axis and use the projection operator (I+Ci,) to project out basis states 
that transform as one of its components. This is equivalent to projecting out the states which 
transform as a given row of the irreducible representation; the latter are not listed in standard 
group theoretical character tables. 

Computation of static properties such as spin densities and spin- spin correlation functions 
in the eigenstates of the Hamiltonian is rendered simple due to the orthogonality of the constant 
Ms basis. The site spin operators such as the z-component of the spin are diagonal in this basis, 
while other operators such as the raising and lowering operators, though not diagonal in this 
basis, have a very simple matrix representation. In computing various properties, the procedure 
we follow is to express the eigenstates in the unsymmetrized constant Ms basis and to compute 
the desired properties using this representation. 

The additional steps involved in the hybrid VB-Constant Ms method are (i) constructing the 
C matrix, whose i th row contains the coefficients of the constant M5 functions appearing in the 
i th VB basis function. This is a very fast step as the constant M5 states are an ordered sequence 
of integers and a VB state with n lines is a linear combination of 2 n constant M5 functions, 
(ii) In the hybrid approach, computation of the C matrix involves the matrix multiplication, 
C(Lr %f r (R)RM) = CR^. The number of arithmetic operations involved is however very small, 
since both C and R^ are sparse matrices. In both constant M5 and hybrid approaches one has 
to obtain the projection operator Pr by retaining only the orthogonal rows of the matrix R^ or 
Qr respectively. Since the number of orthogonal rows in Qr is far fewer than in R^, this step 
is faster in the hybrid approach than in the constant Ms approach by a factor D(Ts)/D(Tm s ), 
where D(Ts) is the dimensionality of the space of the irreducible representation T with spin 
S and D(Tm s ) is similarly the dimension of the space T with constant M5. This advantage is 
largely off- set by the fact that the R^ matrix in constant Ms basis is more sparse than the Qr 




Figure 3: Schematic of Cu$Fe% cluster. Filled and open circles correspond to Fe and Cu (both 
spin- 1/2) sites respectively. Lines represent the exchange coupling between various spin sites. 



matrix in the hybrid approach. Computation of the eigenvalues in the constant Ms approach is 
slower than in the hybrid approach, since D(r M5 )>D(rs) for most S. This advantage may be 
partially off-set because the Hamiltonian matrix in the hybrid approach is usually more dense. 
The memory required for the hybrid approach is not very different from that of constant Ms 
approach since the matrices though smaller in the hybrid approach, are slightly denser. The 
only additional array required in the hybrid approach is the storage of C matrix, which is very 
sparse. The major advantage of the hybrid approach is that we can obtain a far richer spectrum, 
since we are targeting each spin sector separately, unlike in the constant Ms approach. Thus, 
if we can obtain, say 10 states in each S sector of the In spin- 1/2 problem, we would have 
10(rc+ 1) unique states compared to the constant Ms technique where many of these spin states 
would be repeated in different Ms sectors. 



4 Application to CueFe% 

We have applied the above method to model the susceptibility behavior of the molecule [(Tp)$ 
{H 2 0)f 3 Cu I lFe% I (CN) 2 4\ (C/0 4 ) 4 • \2H 2 -2Et 2 ESQ, where Tp stands for hydrotris (pyra- 
zolyl) borate (Fig. 3). In this molecule the Cu 11 ions as well as the Fe 111 are in spin- 1/2 state. 
The eight Fe m ions are at cube corners and the six Cu 11 ions are on the outward perpendicular 
to the face centers of the cube. Each Cu 11 ion is connected to the four nearest Fe 111 ions via 
ferromagnetic exchange interactions. There are no Fe-Fe or Cu-Cu interactions. The molecule 
has a spin 7 ground state. This system has a very high symmetry of the cube and incorporates 
all the complexities that can be encountered in the application of our technique. From the sus- 
ceptibility data, the strength of the exchange interaction J, was estimated to be 30 cm -1 [|27l . 
Because of the high symmetry of this system, we chose this as an example for applying our 



Table 1: Dimensionalities of the total spin spaces of a system of 14 spin-1/2 objects. D(S) is 
the dimensionality of the constant S basis and D(M$) is the dimensionality of the constant Ms 
basis. 
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Table 2: Dimensionalities of various subspaces of the Cu^Fe^ 1 cluster for irreducible represen- 
tation r. 
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technique. The dimensions of the various subspaces are given in Table 1. 

Using the hybrid VB-constant Ms method, we have broken down the space in each total 
spin sector into basis states that transform as different irreducible representations of the cubic 
point group. The dimensionalities of the various subspaces are shown in Table 2. The subspaces 
transforming as the E representations are broken down into subspaces of half their dimension 
by quantizing the system along a C 2 axis. Similarly, the subspaces transforming as the T repre- 
sentations are broken down into one third their dimensions in Table 2, by using a C3 axis as the 
axis of quantization. 

We have set up the Hamiltonian matrix in each of the subspaces and obtained all the eigen- 
values. We have also used a constant Ms basis and using the full cubic symmetry, factored the 
space into various irreducible representations and obtained all the eigenvalues in each subspace. 
From the eigenvectors, we have computed the total spin of the state. We find a one to one cor- 
respondence to numerical accuracy, between the two sets of calculations. We have also fitted 
the %T vs. T experimental plot by using the full spectrum of the Heisenberg Hamiltonian and 
computing ll28l 



10 



o 

E 

=3 

E 

CD 



E 




300 



Figure 4: Fit of the %T vs. T plot for the Cu^Fe^ 1 cluster. The best fit parameters are, J = 27.2 
cm -1 (ferromagnetic), zJ' = -0.008cm _1 (antiferromagnetic), D = -0.15cm _1 and g = 2.1. 
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g 2 F(J,T) 



l-zJ'F(J,T)/k B T 



(9) 



where, we have taken the g factor to be 2.1, the ferromagnetic exchange constant J to be 27.2 
cm -1 . Here, jT is in units of Afy/g. The function F(J,T) is given by, 



F{J,T) 



LsLM s M 2 s expl-Eo(S,M s )/k B T} 
LslM«exp[-Eo(S t Ms)/k B T] 



(10) 



with Eq(S,Ms) being the eigenvalue of the sum of exchange Hamiltonian and the magnetic 
anisotropy term DSf and zJ is the intermolecular exchange interaction. Here we have assumed 
that the molecular anisotropy is along the global z-axis, and this term is treated as a perturbation 
to the exchange Hamiltonian in Eq. 1. In Fig. 4, we show the fit of the experimental data to the 
model. 
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5 Summary and Outlook 

The problem of exploiting total spin invariance together with spatial symmetries, especially of 
non-Abelian point groups has been a long standing one. While full spin symmetry adaptation 
can be achieved by various techniques such as the use of permutation groups, unitary groups 
and the VB method, the last mentioned is the easiest and provides easy chemical visualiza- 
tion of the basis states. The main objection to the VB technique had been that the basis is 
nonorthogonal and leads to nonsymmetric representation of the Hamiltonian matrix. However, 
with the modification of the Davidson's algorithm [|32ll to nonsymmetric matrices by Rettrup 
11331 . this objection has ceased to be important. The ease with which VB states with any given 
total spin can be generated from objects with assorted individual spins, is an advantage which 
far outweighs the other historical objections to the VB method. 

However, even the VB basis suffers from the serious disadvantage like all other spin adapted 
methods, when the question of full spatial symmetry adaptation comes up. The constant A/5 
basis methods do not suffer from this disadvantage. In this paper, we have demonstrated, how 
by combining the ease of spin symmetry adaptation of the VB method with the spatial symmetry 
exploitation of the constant Ms methods, we can devise a scheme which is fully spin and spatial 
symmetry adapted. This has been possible because of the ease of transformation of the VB basis 
to the constant Ms basis. We have demonstrated the power of the method by applying it to the 
exchange Hamiltonian of the molecular magnet C^Fes which has cubic symmetry. We note 
that the hybrid VB-constant Ms method also allows easy manipulation of the eigenstates of the 
Hamiltonian for computing properties. The method described here can easily be extended to 
fermionic systems and should provide a significant improvement for obtaining exact eigenstates 
of spin conserving model Hamiltonians. 
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